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ABSTRACT 

The intergalactic magnetic field (IGMF) may leave an imprint on the anisotropy prop- 
erties of the extragalactic gamma-ray background, through its effect on electromag- 
netic cascades triggered by interactions between very high energy photons and the 
extragalactic background light. A strong IGMF will deflect secondary particles pro- 
duced in these cascades and will thus tend to isotropize lower energy cascade photons, 
thus inducing a modulation in the anisotropy energy spectrum of the gamma-ray back- 
ground. Here we present a simple, proof-of-concept calculation of the magnitude of 
this effect and demonstrate that the two extreme cases (zero IGMF and IGMF strong 
enough to completely isotropize cascade photons) would be separable by ten years of 
Fermi observations and reasonable model parameters for the gamma-ray background. 
The anisotropy energy spectrum of the Fermi gamma-ray background could thus be 
used as a probe of the IGMF strength. 
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1 INTRODUCTION 

Cosmic magnetic fields are expected to play a fundamental 

role in the physics of a large variety of astrophysical sys- 

tems. The fate of high energy particles, such as ultra-high 

energy cosmic rays, as they propagate through the Universe 

and the acceleration of charged particles in astrophysical 

objects hinge, in large part, on magnetic fields within their 
sources and those permeating the Universe. Large-scale mag- 
netic fields such as those found in galaxies and clusters of 
galaxies have been a mystery for the past several decades. 
While they are generally thought to be the result of the 

amplification of weak seed fields, the origins of these seed 
fields, whether they be cosmological or astrophysical in na- 
ture, are largely unknown. A definitive measurement of the 

intergalactic magnetic field (IGMF) could provide a funda- 
mental step in resolving the questions of the origins of cos- 
mic magnetic fields and their impact on the evolution of 
the systems in which they reside, but sufficiently constrain- 
ing observations have thus far remained elusive. Previously, 
upper bounds of ~ 10~~ 9 G have been found through Fara- 
day rotation limits of polarized radio emission from distant 
quasars and the study of the effect of magnetic fields on 
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the anisotropy of the cosmic microwave background (for a 
brief review of observational and theoretical techniques to 
constrain the IGMF, see e.g., Neronov & Semikoz 2009, and 
references therein). The recent availability of data from the 
Fermi Large Area Telescope (LAT) has renewed interest in 
the IGMF as observations of the gamma-ray sky and its par- 
ticipating source populations provide a unique perspective 
on the strength and structure of the IGMF. 

The gamma-ray sky consists of resolved point sources 
(such as normal galaxies and active galaxies, pulsars, etc.), 
transient gamma- ray sources (e.g., gamma-ray bursts), and 
the diffuse gamma-ray radiation comprised of emission from 
the Galaxy and the isotropic (presumably, extragalactic) 
gamma-rav background (EGB). The EGB is a window into 
the high-energy processes in the Universe the origins of 
which have been the subject of much debate; however, it is 
expected that emission arising from unresolved, extragalac- 
tic point sources, such as blazars and star-forming galax- 
ies, comprises a sizable contribution to the EGB (see e.g., 
Stecker & Salamon 1996; Fields et al. 2010: Abdo et. al. 
2010d; Stecker & Venters 2011, and references therein). Ad- 
ditionally, many of these extragalactic point sources are also 
sources of very-high energy (VHE) gamma rays, 1 which in- 


* In this paper, we take VHE to be ~ TeV. 
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teract with the soft photons of the extragalactic background 
light (EBL), consisting of the cosmic microwave background 
(CMB) and the background of infrared, optical, and ultra- 
violet radiation from direct and dust-reprocessed starlight. 
The interactions of VHE photons with the EBL initiate 
electromagnetic (EM) cascades, giving rise to another con- 
tribution to the EGB (see e.g., Coppi & Aharonian 1997; 
Kneiske & Mannheim 2008; Inoue & Totani 2009; Venters 

2010) . However, the degree to which astrophysical sources 
contribute to the EGB depends on their distribution with re- 
spect to luminosity and redshift (the gamma-ray luminosity 
function, GLF), their intrinsic spectra at gamma-ray ener- 
gies, and the nature of the EBL, all of which remain the sub- 
jects of intense debate (see e.g., Venters 2010; Fields et al. 
2010; Abdo et al. 2010b, d; Makiya et al. 2011; Inoue 2011; 
Stecker & Venters 2011; Abazajian et al. 2011). 

A promising complimentary technique for investigating 
the contributions of astrophysical sources to the EGB in- 
volves studying its intensity anisotropy properties as a func- 
tion of energy. As pointed out in Siegal-Gaskins & Pavli- 
dou (2009; see also Hensley et al. 2010; Siegal-Gaskins et al. 

2011) , the contribution of a given class of sources to the 
overall fluctuation angular power at a given angular scale of 
the EGB at a given energy is weighted by the square of its 
fractional contribution to the intensity of the EGB at that 
energy. As such, if the relative contributions of the various 
EGB contributors change with energy, so too will the an- 
gular power as a function of energy. Expectations for the 
anisotropies of the commonly invoked astrophysical contrib- 
utors to the EGB are that star-forming galaxies contribute 
very little anisotropy owing to the fact that they are so nu- 
merous and individually faint, while blazars, being far more 
rare than star-forming galaxies but much brighter, could 
have a large contribution to the EGB anisotropy (Ando et al. 
2007; Ando & Pavlidou 2009; Siegal-Gaskins & Pavlidou 
2009; Hensley et al. 2010; Cuoco et al. 2011). It has also 
been suggested that annihilation in dark matter substruc- 
ture in the Galactic halo could also have a large contribu- 
tion to the EGB anisotropy, even if the contribution to the 
EGB intensity from dark matter annihilation is subdom- 
inant with respect to other contributions (see e.g., Ando 
et al. 2007; Siegal-Gaskins 2008; Taoso et al. 2009; Ando 
2009; Siegal-Gaskins & Pavlidou 2009: Hensley et al. 2010). 
Early anisotropy results from the Fermf-LAT Collaboration 
(Cuoco et al. 2011) suggest that the energy dependence of 
the intensity anisotropy of the EGB may be consistent with 
arising from a source population with spectra similar to that 
of blazars (a ~ 2.4), even though blazars may comprise 
< 30% of the EGB intensity (Abdo et al. 2010d). Thus, 
even contributions to the EGB that are subdominant with 
respect to other contributions can have a profound impact 
on the anisotropy of the EGB, especially if its fractional 
contribution changes with energy. As such, it could prove a 
powerful tool in detecting changes in the EGB anisotropy 
resulting from the magnetic deflection of EM cascades from 
blazars. 

In interacting with EBL photons, VHE gamma rays 
initiate EM cascades of photons, electrons, and positrons. 
The interaction of a VHE gamma ray* with an EBL photon 
produces a pair of electrons and positrons, which will, in 
turn, Inverse Compton scatter EBL photons to high ener- 
gies, These up-scattered photons will, in turn, pair produce, 


and the process continues until the energies of the resulting 
photons are low enough that pair production is no longer ef- 
ficient. In the presence of an IGMF, the charged particles in 
the cascade are deflected away from the direction of propaga- 
tion of the primary photon. As these charged particles prop- 
agate through the EBL, they can up-scatter soft photons 
into the line of sight of the observer, though they will appear 
to be coming from a different direction. Observationally, 
this process results in a “halo” of gamma rays around the 
source. The detection of a halo around a gamma-ray point 
source would not only provide strong evidence for an IGMF, 
but would also provide an indication of the strength of the 
IGMF. Several groups have already conducted searches for 
extended emission around gamma-ray sources (see e.g., Aha- 
ronian et al. 2001; Neronov & Vovk 2010; Ando & Kusenko 
2010; Aleksic et al. 2010; Neronov et al. 2011), though as yet, 
no significant indication of gamma-ray halos distinguishable 
from instrumental effects have been reported. However, the 
absence of significant halo emission could imply a strong 
IGMF that would extend the EM cascade halo such that the 
surface brightness of the extended source falls below the de- 
tection threshold (Neronov & Vovk 2010; Ando & Kusenko 
2010). In addition to possibly creating a halo, the deflection 
of EM cascades by the IGMF could result in a reduction 
in observable cascade radiation. Analyses employing simul- 
taneous observations of blazars by Atmospheric Cherenkov 
Telescopes and Fermi make use of this fact to derive limits 
on the IGMF (Dermer et al. 2011; Essey et al. 2011; Tay- 
lor et al. 2011). These techniques are quite promising and, 
with more data from Fermi , will continue to shed light on 
the IGMF. In this paper, we consider a different promising 
technique that makes use of the anisotropy energy spectrum 
of the EGB to constrain the IGMF. 

For any cosmological population, such as blazars, that 
emits gamma rays at very high energies, the effect of EM 
cascading results in a flux suppression at the highest en- 
ergies and enhancement at lower energies as the cascades 
redistribute radiation from the higher energies to the lower 
energies (Coppi & Aharonian 1997; Kneiske & Mannheim 
2008; Inoue & Totani 2009; Venters 2010). In considering 
the spectral subpopulations of blazars, flat-spectrum radio 
quasars (FSRQs) and BL-Lacertae type objects (BL Lacs), 
the intrinsic spectra of FSRQs is much softer than that of 
both BL Lacs and EM cascades, and their relative promi- 
nence with respect to the EGB can change with energy (Ven- 
ters 2010; Venters & Pavlidou 2011). As such, the fluctua- 
tion angular power in the EGB should change with energy 
even if FSRQs and BL Lacs have similar anisotropy prop- 
erties. Furthermore, as noted earlier, the deflection of EM 
cascades by the IGMF can spatially extend the emission 
around gamma-ray sources and reduce the observable cas- 
cade radiation. Thus, the impact of a strong IGMF is to 
reduce the blazar anisotropy both because of the increased 
angular scale of the sources and because of the reduction in 
the overall blazar contribution to the EGB (assuming that 
blazars are the most anisotropic EGB component). On the 
other hand, in the absence of a strong IGMF, the EM cas- 
cades are highly collimated in the direction of the primary 
photon (5 ~ m e /E e ; Neronov & Semikoz 2009), and their 
impact on the blazar anisotropy is to reinforce it due to the 
enhancement in the overall blazar contribution to the EGB. 
Hence, in the two extreme cases for the IGMF, EM cascades 
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could have a profound and opposite impact on the blazar 
anisotropy and, by extension, the EGB anisotropy. There- 
fore, studies of the EGB anisotropy could prove a sensitive 
probe into the nature of the IGMF. 

In this paper, we explore the impact of EM cascade ra- 
diation from blazars on the IGRB anisotropy in the extreme 
cases of a negligible IGMF and a strong enough IGMF to 
essentially isotropize the cascade emission (i.e., from the ob- 
servationally practical standpoint). This study is meant to 
serve as a proof of concept for the feasibility of using IGRB 
anisotropies as an IGMF probe. In Section 2, we present 
the methodology and parameters employed to determine the 
anisotropy energy spectrum and the EM cascades. In Sec- 
tion 3, we present the results of the calculation. In Section 
4, we discuss our results and conclude. 


2 METHODOLOGY & INPUTS 

2.1 Collective Intensity of Unresolved 
Astrophysical Sources 


The contribution to the EGB due to unresolved blazars can 
be viewed as the superposition of the collective intensity of 
intrinsic blazar spectra and the contribution from the cas- 
cade radiation from the interactions of VHE photons from 
blazars with the EBL: 

lE(Eo)=I i B (Eo) + I%(E 0 ), (1) 


where the intensity, Ie(Eq), is given in units of photons 
per unit area per unit time per unit solid angle per unit 
energy emitted at observer frame energy, Eq. In calculating 
the collective spectrum of unresolved blazars, we follow the 
formalism as outlined in Venters & Pavlidou (2011). Blazar 
spectra are taken to be smoothly broken power laws: 


F e (Eo) = F 0 




(2) 


where Fe(Eq) is the differential photon flux in units of pho- 
tons per unit area per unit energy per unit time, Eb is 
the break energy, a\ = a — Aon is the low-energy slope, 
ftj = a + Act 2 is the high-energy slope, a is the gamma-ray 
photon spectral index of the blazar assuming an unbroken 
power- law spectrum, 2 and n quantifies the sharpness of the 
transition from the low-energy power law to the high-energy 
power law (taken to be 1). 

The total flux, F, of photons with energies greater than 
some fiducial energy, Ef (taken to be 100 MeV), is found by 
integrating Fe(Eo) over energy, 


F — Fo 




( 3 ) 


Then, the contribution of a single unresolved blazar to the 
EGB is 


F UEo/Eb)^ + (E 0 /E h )° 


l-r./^ l(Eo/E b )° 


(E 0 /E b )^r dE 0 


( 4 ) 


where the flux of one source is uniformly distributed over 


2 For the sake of clarity, we use a for the photon index rather 
than F, which is commonly used in the literature. 


the entire sky in anticipation of an isotropically distributed 
cosmological population. Including absorption due to the 
EBL, the contribution becomes 

T = F[(E 0 /E b Y' + (E 0 /E b rT l -rtish.z) (5) 
1 4n [(E 0 /Eb) ai + (E 0 /E b ) a 2] _1 dE 0 1 

where t(Eo,z ) is the optical depth for a photon with ob- 
server frame energy, Eo, emitted at redshift, z. In determin- 
ing the total unresolved blazar contribution including EBL 
absorption, one would generally rewrite Equation 5 in terms 
of the blazar gamma-ray luminosity and redshift, multiply 
by the blazar GLF and intrinsic spectral index distribution 
(corrected for the spectral bias endemic to flux-limited sur- 
veys 3 as per Venters et al. (2009)), and finally integrate over 
gamma-ray luminosity, redshift, and spectral index (see e.g., 
Venters et al. 2009). However, since the focus of this paper 
is to study the relative changes with energy in the EGB 
anisotropy, we do not fully calculate the unresolved blazar 
contribution to the EGB. Rather, we assume typical lumi- 
nosities and redshifts for all blazars. As such, the only inte- 
gral that must be performed is that over the spectral index: 


Ie(Eo) 


lol 




S(E f ,a) 


(6) 


where I o is a normalization constant, p(a) is the measured 
spectral index distribution (SID), and 



( 7 ) 


We fix the normalization constant such that the total astro- 
physical (blazars, star-forming galaxies, and EM cascades) 
contribution to the EGB fits the spectrum of the EGB as 
measured by the Fermi- LAT (Abdo et al. 2010c). In deter- 
mining the measured SID, we make use of the results pre- 
sented in Venters & Pavlidou (2011) from a likelihood analy- 
sis fitting Fermi-LAT First Catalog (Abdo et al. 2010e) FS- 
RQs and BL Lacs to Gaussian SIDs accounting for errors in 
measurement of individual blazar spectral indices. For these 
blazar subpopulations, the maximum-likelihood Gaussian 
SIDs can be characterized by means (ao) and spreads (err,) 
with parameters determined to be ao = 2.45 and ao = 0.16 
for FSRQs and ao = 2.17 and ao = 0.23 for BL Lacs. Fol- 
lowing Venters & Pavlidou (2011), we model the spectral 
breaks by taking Aon = 0.1, A «2 = 0.9, and E b ,o = 4 GeV 
for FSRQs. For BL Lacs, we take Eb, o = 10 TeV. Since the 
AGN population peaks ~ 1, we take the typical redshift of 
blazars to be z src — 1. In modeling the suppression of the 
source brightness due to magnetic deflection of cascades (see 
Appendix), we take the jet opening angle of blazars to be 
~ 1°. The model of the EBL was taken from Franceschini 
et al. (2008). 

We also include a model of the star-forming galaxy (SF) 
contribution to the EGB. In this case, we assume that the 


3 We note that the Fermi survey is not exactly a flux-limited 
survey due to the non-uniformity of the total diffuse background 
throughout the sky. However, for the purposes of this paper, we 
neglect this effect. 
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gamma-ray luminosity of an SF galaxy is mainly due to 
the decay of neutral pions created through galactic cosmic- 
ray interactions with neutral hydrogen in the galaxy. It is 
thought that galactic cosmic rays are accelerated in su- 
pernova remnants and then diffuse throughout the galaxy. 
Hence, the flux of galactic cosmic rays is expected to be 
proportional to the galaxy’s supernova rate, which is, in 
turn, proportional to its star-formation rate. Furthermore, 
the neutral hydrogen gas provides the fuel for forming stars. 
As such, the expectation is that an SF galaxy’s gamma- 
ray luminosity can be parameterized in terms of its star- 
formation rate (Pavlidou & Fields 2002; Abdo et al. 2010a; 
Fields et al. 2010; Makiya et al. 2011; Stecker & Venters 
2011). For the purposes of this paper, we model the star- 
forming galaxy contribution using the infrared luminosity 
function model of Stecker & Venters (2011). In this model, 
the star-formation rate of an SF galaxy is assumed to be 
proportional to its infrared luminosity; thus, the gamma- 
ray luminosity of an SF galaxy is proportional to a power 


of its infrared luminosity. In Stecker & Venters (2011), the 
relationship between the gamma-ray luminosity of an SF 
galaxy and its infrared luminosity was determined by fitting 
the measured gamma-ray luminosities of Fermi-detected SF 
galaxies to their infrared luminosities. The SF galaxy con- 
tribution could then be determined by integrating over an 
infrared luminosity function of SF galaxies, which they took 
from Hopkins et al. (2010). 

2.2 Electromagnetic Cascades from Blazars 

The cascade intensity is given by 


E dtdAdQ.dE ’ 

where dN c /dE is the spectrum of cascade photons due to 
pair production and Inverse Compton scattering (for full 
equations and derivations, see Venters 2010): 


dN c 

dEn 


(So) 


T 


(l+z) 


drdSLy 
dzdE , 


P(f', E p ,z) 


dNr-t (Fo(l + z)) 


dE 


dNr 2 (E 0 (l + z )) 
d.E 


-t(Eq.z) 


dE p dz 


(9) 


where dNr/dE is the spectrum of Inverse Compton scat- 
tered radiation per electron of Lorentz factor, F, P(f\ E p , z) 
is the probability that the pair production interaction at a 
given redshift, z, of a primary photon of energy E p will pro- 
duce electron-type particles of energies E e i = / x E p and 
E e 2 = (l-f)xEp, Fi = fxEp/mc 2 , T 2 = (1 -f)xE p /m.c 2 , 
and d 2 Nj /dzdEp is the continuous spectrum of photons un- 
dergoing pair production interactions as a function of red- 
shift and primary energy. 

The full numerical integration of Equation 9 is compu- 
tationally cumbersome. As such, in order to calculate the 


cascade spectrum, one typically makes use of a Monte Carlo 
code, such as Cascata (Venters 2010), that propagates pho- 
tons from known spectra and redshift distributions and cal- 
culates the cascade spectrum of each photon. However, for 
the purposes of this paper, we employ a more semi-analytical 
model of the cascade development. Since we have assumed a 
typical redshift for all blazars, we have eliminated the need 
for randomly generated source redshifts. Instead, the spec- 
trum of pair-producing photons, d 2 N~/dzdE p , can readily 
be determined from the absorbed collective intensity of un- 
resolved blazars binned in energy: 



l£ s (E 0 )dE 0 = Iq 



dav(a) liEo/Eb)ai+(Eo/E ^ 2] 
P[ ’ S(E/,a) 


(l - e - T(Bo ’ z) ) dE 0 . 


(10) 


For each energy bin, we determined a characteristic red- 
shift for the start of the cascades based on the criterion 
that the characteristic redshift is that at which the num- 
ber of photons that have been “absorbed” in propagating 
from the source redshift is half the number of the total pho- 
tons “absorbed” in propagating to the present epoch. Using 
this criterion, the characteristic redshift can be calculated 
numerically from 


r,(Ea,z r ) — In 


[1 

i 2 


^ e ’H- E 0,Zsrc) 



(ID 


Also, we assume that each secondary electron produced in a 
pair-production interaction carries half of the energy of the 
primary photon. Finally, since the number density of CMB 
photons is much higher than those of the other wavelengths 
that comprise the EBL, we assume that the seed photons for 


the Inverse Compton interactions are CMB photons. In so 
doing, we can neglect Klein-Nishina effects in determining 
the spectrum of scattered radiation. The spectrum of scat- 
tered radiation per electron of Lorentz factor, F is given by 
(Venters 2010) 


dN r _ 

f d3Nr ' dtde- 

f d 3 N r ' 

dt 

de 1 

f J j ■ (vvwC — — 

/ dtdeide 

1 dtdeide 

dE e 


(12) 


where F' is the electron Lorentz factor at time, t, e is the 
initial soft photon energy, ei is the scattered photon energy, 
d 3 N r > /dtdeide is the differential scattered photon spectrum 
given by (Blumenthal & Gould 1970) 


d 3 N r , 

dtdeide 


3(Tt c n(e) 
16F 4 


2ei In 


ei 

4F' 2 e 


• ei • 


2r ,2 e 


(13) 
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n(e)de is the number density of soft photons, and j dE e /dt 
is the rate of energy loss of the electron given by 


d.E e 

dt 


4 tV2 
CT'1'Cl 

3 


£ 


i(e)de . 


(14) 


We note that this formalism assumes continuous energy 
losses for the electrons, which is appropriate for the regime 
that we are considering. Furthermore, this simplified ap- 
proach has the advantage that recycling cascades (in which 
Inverse Compton photons up-scattered to high energies are 
also allowed to initiate EM cascades) is more feasible than 
in the full Monte Carlo approach. As such, for BL Lacs, 
we also calculated the spectrum of second generation cas- 
cades. Notably, the effect of EM cascades is to redistribute 
high-energy radiation to lower energies; as such, the high- 
energy component of second generation cascades necessary 
to produce higher generation cascades is greatly reduced. 
Thus, the contribution from generations beyond the sec- 
ond generation is negligible. For FSEQs, the intrinsic spec- 
tra of the sources are sufficiently soft that even the first 
generation of cascades makes a small contribution to their 
overall collective intensity. 4 The suppression of the source 
brightness due to magnetic deflection of cascades (see Ap- 
pendix) is also calculated, and we assume a uniform mag- 
netic field that evolves solely due to the expansion of the 
Universe [ B(z ) ~ Bq(\ + z) 2 ; Neronov & Semikoz 2009], 
where Bo ~ 5 x 10~ 14 G. We should also note that for the 
purposes of this paper, we assume that the correlation length 
of the magnetic field is much larger than electron cooling 
distances. Alternative scenarios for the magnetic field and 
source spectra will be explored in future publications. 


2.3 Anisotropy Energy Spectrum 

The anisotropy energy spectrum, CpE), is defined as the 
intensity fluctuation angular power (in units of sr) at a given 
angular scale, l, as a function of energy. In the simple case 
of a two-component background, the total angular power is 
given by (Siegal-Gaskins & Pavlidou 2009) 

Cj°* = fic f> + /|cf 5 + 2/ 1 / 2 C'i lx2) , (15) 

where C ( e n ' is the angular power spectrum of component 
(n), f n = In(E)/I to t(E) is the energy-dependent fraction of 
the total emission arising from component ( n ), and C] nxm) 
is the cross-correlation term for components (n) and (rn). 
As is a measure of angular fluctuations in units of the 
mean, it is independent of energy for a single population 
of sources with identical spectra. If the two components are 
uncorrelated, c£ nXm) = 0, and Equation (15) is further sim- 
plified, 

CT = fiC™ + /fcj 2) . (16) 

If additional components of negligible angular power con- 
tribute to the total background signal, then these will not- 
result in additional terms in Equation (16), but rather they 


4 We note that the small amount of cascade radiation is not sim- 

ply the result of the low break energy chosen for this analysis: 
Venters (2010) modeled FSRQ spectra as unbroken power laws 
up to ~ 100 TeV and also found that the FSRQ cascade contri- 
bution was small. 


will affect C\ ot (E) by changing the fractional contributions, 
/, of each of the components that do contribute to the an- 
gular power. 

In the simple scenarios we explore, we always consider 
two components that do contribute to the total angular 
power: BL Lacs and FSRQs. Because in our adopted model 
the two contributions do not add up to the total EGB signal 
measured by Fermi (Abdo et al. 2010c), we assume that the 
rest of the EGB signal originates from a very-low Ce popu- 
lation, such as star-forming galaxies. For simplicity, we take 
C£ F ' = 0 for this latter component. As discussed in the 
previous sections, the energy spectrum of the primary emis- 
sion from FSEQs is very soft at very high energies, so the 
cascade signal from FSEQs is negligible; all cascade emis- 
sion that we consider comes from reprocessed BL Lac VHE 
photons. 

In the limit of low B, there is little angular spread in 
the cascade photons, so these photons appear to originate 
from the same sources as the primary emission. In this case, 
the only effect of the cascades is to alter the energy spec- 
trum of BL Lacs, leaving their anisotropy properties un- 
changed. For this reason, in this scenario, we add the cas- 
cade radiation from BL Lacs to their primary emission in 
f '2 of Equation (16) while considering only primary emis- 
sion from FSRQs in /i; hence, the fractional components 
of Equation (16) become fi(E) = Ifsrq(E) / Itot(E) and 
MB!) = [Ibllzc(E) + /ease (£)] //tot (£), while the cf SRQ) 
and (M BljLac ! remain unchanged. 

In the limit of strong B (Bo ~ few x 10"~ 14 G), we deter- 
mined the scattering angles of the cascades (see Appendix) 
and found that for observed photon energies less than ~ 100 
GeV, the scattering angles were sufficiently large such that 
within the observed angular extent of the cascade emission, 
the number of sources above the Fermi- LAT sensitivity is 
sufficiently large to introduce source confusion 0 . Thus, we 
take the cascade signal as having angular power, C'£ asc) = 0. 
As such, in this scenario, we simply add the cascade emission 
to /tot, while the numerators of fi and / 2 in Equation (16) 
still correspond to the primary (appropriately attenuated) 
emission from FSRQs and BL Lacs, respectively. 

We have considered a range of angular scales between 
l = 50 and l = 250 (l = 150 with A£ — 100). For simplicity, 
we have taken Cf SRQ) = C< BLLac) = 6 x 10~ 5 sr, which 
in both scenarios, gives C\ ot ~ 10~° sr around 1 GeV, con- 
sistent with the Fermi measurement of angular anisotropies 
in the EGB at these energies (Cuoeo et al. 2011). We have 
calculated error bars using (Knox 1995) 

= / (2/+ 1)A// Sky ( Ce + Wf ) ’ (17) 

where Cn = (Anfsky/Ny) is the photon noise angular power, 
and We = exp(— Perl/ 2) is the window function of a Gaus- 
sian beam. We have taken f s ^ y = 0.32 (the unmasked frac- 
tion of the sky that would presumably be used for the 
anisotropy data analysis), <?b — o r 6s%/l-51 (where a&%% is the 
68% containment angle of the Fermi-LAT reconstruction), 
and we calculate N~, = dE~A c s{E) fH 0 bs from the 


5 Since the blazar contribution to the EGB is the collective emis- 
sion of unresolved blazars, this criterion is quite conservative. 
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(a) 


(b) 


Figure 1. a: The spectrum of the modeled EGB (thick solid line) together with individual components: intrinsic blazar emission, 
including both FSRQs and BL Lacs (thin solid line); SF galaxies (dashed line); electromagnetic cascades for zero IGMF (dot-dashed 
line) and nonzero IGMF (double dot-dashed). For reference, the spectra of the EGB based on both Fermi (filled circles; Abdo et al. 
2010c) and EGRET (open squares; Sreekumar et al. 1998; open triangles; Strong et al. 2004) are also plotted. b\ The anisotropy energy 
spectrum of the modeled contributions to the EGB. The solid line is the total model assuming zero IGMF. The double dot-dashed line 
is the total model assuming nonzero IGMF. For reference, uncertainties calculated assuming five (dashed error bars) and ten (solid error 
bars) years of observation time are also plotted for four energy bins. 


Fermi-measured EGB intensity (Abdo et al. 2010c), taking 
the Fermi effective area, A e g(E), and a 6S % from the LAT 
performance curves 6 and assuming an observation time, t 0 b s , 
for Fermi of five and ten years. 


3 RESULTS 

In Figure 1(a), we plot the spectra of the modeled contri- 
butions to the EGB. The “total” contribution is the sum 
of the contributions from blazars, SF galaxies, and one of 
the cascade models (the total is similar for both models). 
The blazar component hardens above ~ few GeV due to 
the increased contribution of BL Lacs with respect to FS- 
RQs, as BL Lac spectra remain hard while FSRQs soften. 
The cascade spectrum in the case of the zero IGMF is quite 
flat ( dN/dE oc E~ 2 ) demonstrating the cascade effect of 
redistributing radiation from high energies to lower ener- 
gies. The sum of the three components fits quite well the 
spectrum of the EGB as measured by Fermi. Notably, the 
spectra of the cascades for the two cases of the IGMF while 
similar at the highest energies are quite dissimilar at the 
lower energies. This is due to the loss of some cascade ra- 
diation in the strong IGMF case as significant amounts of 
the cascades are deflected away from the line-of-sight of the 


6 http : //www . slac . Stanford . edu/exp/glast/groups/canda/ 
lat_Perf onaance . htm 


observer. However, it should be noted that given our sim- 
ple model of the blazar contribution to the the EGB, we 
were only able to account for cascade losses due to deflec- 
tion. That is, in our simple picture, we are calculating the 
amount of radiation that an observer looking head-on at a 
given source would be able to detect. It is likely that just as 
cascades from one source are deflected out of the observer’s 
line-of-sight, cascades from other sources that would oth- 
erwise be undetectable (e.g., misaligned jets) are deflected 
into the observer’s line of sight. As such, the cascade inten- 
sity plotted for the strong IGMF case should be viewed as 
a lower limit. Nevertheless, we note that even if the strong 
IGMF cascade radiation were comparable to that of the zero 
IGMF case, the anisotropy would still be different since the 
angular properties of the sources would be different. 

In Figure 1(b), we plot the anisotropy energy spectra of 
the EGB for the two cases of the IGMF. Notably, there is an 
appreciable difference between the two cases starting at ~ 
few GeV. This is due to the fact that at these energies, the 
contribution from EM cascades relative to the other compo- 
nents is much more significant than at lower energies. As a 
result, the isotropization of the EGB due to cascades in the 
case of the strong IGMF is much more significant. Aiding in 
the distinction is that the cascades augment the anisotropy 
in the case of zero IGMF, because they enhance the contri- 
bution of the most-anisotropic component. Also plotted are 
uncertainties for four bins in energy calculated assuming five 
and ten years of observation time. The uncertainties indicate 
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that perhaps within five, and certainly within ten, years of 
observations, Fermi will be able to distinguish between at 
least these two extreme cases. 


4 DISCUSSION & CONCLUSION 

We have studied the effect of a strong IGMF on EM cas- 
cades from blazars and the resulting impact on the cascade 
component of the EGB and the anisotropy energy spectrum 
of the EGB. We have shown that while the spectrum of cas- 
cade radiation in the case of zero IGMF is flat in E 2 , the 
spectrum in the case of a strong IGMF is suppressed at low 
energies, though there could be some additional contribu- 
tion arising from radiation deflected into the line-of-sight 
of the observer. Interestingly, in such a case, radiation that 
was initially emitted away from the observer and would not 
be observable in the case of zero IGMF could, in princi- 
ple, be observable in the nonzero IGMF scenario. As such, 
a definitive indication of a nonzero IGMF could have pro- 
found implications for the components of the EGB. We will 
address this possibility in a future publication. 

For the two extreme cases of IGMF strength we have 
considered here, we have calculated the anisotropy energy 
spectrum of the EGB. In the case of zero IGMF, the ef- 
fect of cascades is to augment the anisotropy arising from 
blazars since the cascades increase the contribution to the 
EGB arising from blazars as they are highly collimated with 
the primary emission. In the case of the strong IGMF, the 
cascades are deflected substantially and the angular proper- 
ties of the parent blazar population are blurred. Thus, in the 
case of a strong IGMF, the effect of cascades is to reduce the 
anisotropy of the EGB as they hinder the anisotropy arising 
from blazars rather than reinforcing it. At energies at which 
the relative contribution to the EGB arising from EM cas- 
cades is significant, the anisotropy energy spectra of the two 
cases become distinguishable. We have also calculated the 
uncertainties in the anisotropy energy spectra in four bins 
of energy assuming observation times of five and ten years 
with Fermi. Within a few short years, Fermi will be able to 
use this technique to distinguish between the two cases and 
constrain the IGMF. 

We note that an advantage of our approach is that 
rather than focusing on a few select sources, our method em- 
ploys information from the whole observable sky. As such, 
our approach could prove capable of providing more sensitive 
constraints on the global properties of the IGMF. Further- 
more, our approach is not limited to uniform magnetic fields 
with large coherent lengths; other field configurations can be 
probed with this technique. In practice, in order to use this 
technique to provide an IGMF measurement, intermediate 
cases of IGMF strength need to be considered as well, and 
the uncertainties entering the calculation of the expected 
anisotropy energy spectra in each case carefully evaluated. In 
the case of intermediate IGMF strengths, a numerical sim- 
ulation of the expected angular power is appropriate, as the 
cascade radiation will have somewhat, but not completely, 
suppressed angular anisotropy at small scales compared to 
the primary BL Lac emission, -which will in addition be 
correlated with the latter. As such, for intermediate IGMF 
strengths, a simple analytical calculation of the cascade ef- 
fect on the total angular power is not possible. Concerning 


uncertainties, these enter our calculation primarily through 
the model assumptions for each EGB component (which can 
in turn self-consistently determine their anisotropy proper- 
ties, a calculation which we have not performed here), and 
through the EBL model. We will expand our investigation 
of the practical usage of the anisotropy energy spectrum of 
the EGB to constrain the properties of the IGMF in future 
publications. 

In the meantime, what this simple, proof-of-concept cal- 
culation has already demonstrated is that the anisotropy 
energy spectrum can be a powerful tool in constraining the 
IGMF. With a few more years of data from Fermi, we can be- 
gin to probe the IGMF and uncover more clues into the ori- 
gins of ultra-high energy cosmic rays and large-scale struc- 
ture formation. 
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APPENDIX A: SUPPRESSION OF THE 
SOURCE BRIGHTNESS DUE TO MAGNETIC 
DEFLECTION OF EM CASCADES 


In calculating the suppression of the source brightness due to 
magnetic deflection of EM cascades, we must determine the 
angular spread of the cascades, 9 stc , as viewed by a funda- 
mental observer at the source redshift, z src (Figure Al(a)). 
In order to determine 0 src for multiple cascades, we begin 
with considering the problem for one generation of cascades 
(Figure Al(b)). In this case, 9i, the angular spread of the 
cascade as viewed by a fundamental observer at z\, can be 
expressed in terms of quantities evaluated at the redshift 
of interaction of the “primary” photon, Z2, and the angular 
diameter distance from z src = zi to Z2, D A {z\,Z2): 


d-2 


D a (z 


(Al) 


where di is the proper lateral distance the electron has trav- 
eled as measured by a fundamental observer at z-j (see e.g.. 
Longair 2008), and we are neglecting Klein-Nishina effects, 
so the electron cooling distance, dfj, is not cosmological (i.e., 
za « 22 ). By geometry, do is given by 

di — 1 2 sin d 2 = 2r2 sin* do , (A2) 


where is the straight path of d|, r 2 is the Larmor radius 
at zo of an electron with energy, E%, in magnetic field of 


strength, R(za) = i?o(l+Z 2 ) 2 , given by (Neronov & Semikoz 

2009) 


r 2 = (3 x 10 28 cm) 


B(z 2 ) 
10" 18 G 


El \ 

10 TeV J 


(A3) 


and 62 — (1/2 ) 91 is angle through which the electron is 
deflected away from the straight-line path of the “primary” 
photon, which is equal to half of the magnetic deflection 
angle of the electron given by (Neronov & Semikoz 2009) 



(3 x 10~ 6 ) (I + Z 2 )” 


Bq 


IQ** 18 G 


m v 

10 TeV J 


(A4) 


The angular diameter distance from Zi to Z2 can be found 
from the angular diameter distance from Z2 to Zi through 
the Reciprocity Theorem: 


D A ( Zl ,z 2 ) = 1±I1Da(z 2 ,z 1 ), (A5) 

1 T ^2 

where Da ( 22 , z\) for a flat universe is given by 


Da (z 2 ,zi)= — [D e ( 2 ! ) - D c (z 2 )] , (A6) 

1 + Zl 

and D c (z) is the comoving distance at z. 

Now that we have 9 1 in terms of known quantities, we 
can find 9 src for two generations of cascades using an equa- 
tion similar to Equation Al. Returning to Figure Al(a), 
we note that in this case, the primary photon is emitted 
at Zsrc and interacts at 21 emitting electron-type particles 
that do not propagate cosmological distances. One electron- 
type particle is deflected away from the straight-line path of 
the primary particle through an angle, <5i, and up-scatters 
CMB photons. The up-scattered photons then propagate un- 
til they interact at 22 , creating more electron- type particles. 
Since we are neglecting generations of cascades higher than 
the second generation, we have 9\ from the formalism de- 
veloped above. The angle through which the electron at z\ 
is deflected away from the straight-line path of the primary 
photon, <5i, is also given by Equation A4, but for quantities 
evaluated at z \ . 

Having found 9\ and <5i, we have the total angular 
spread of the cascades as viewed by a fundamental observer 
at zi, 9\ = 9i + 61. The lateral spread of the cascade as 
measured by a fundamental observer at 22 , D2, is given by 
Equation Al: 


D2 = 6\D A { Zl ,z 2 ). 


(AT) 


Finally, we can again use Equation Al to find the angular 
spread of the cascades as viewed by a fundamental observer 
at the source redshift: 


9$rc — 


P2 

D A {z S tc,Z2) ' 


(A8) 


The flux suppression of the cascade due to magnetic deflec- 
tion can be found from the flux of the cascade in the absence 
of magnetic fields by relating the solid angle of the deflected 
emission, A Qd, to that of the unaffected emission (in this 
case, the solid angle of the jet), Ail: 

Af> 

F &cid = ■ (A9) 


We note that this procedure is only valid for non- 
isotropic sources. Conceptually, there should be no flux sup- 
pression for an isotropic source because the symmetry of the 
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Figure 1 . a: Geometry for the angular deflection of two generations of EM cascades, b: Geometry for the angular deflection of one 
generation of EM cascades. 


problem demands that whatever emission that is deflected This paper has been typeset from a TgX/ I^Tl^X file prepared 
out of the observer’s cone of sight is replaced by emission by the author, 
being deflected in. 
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